Parity effect in superconducting aluminum single electron transistors with spatial gap 

profile controlled by film thickness 
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We propose a novel method for suppression of quasiparticle poisoning in Al Coulomb blockade 
devices. The method is based on creation of a proper energy gap profile along the device. In contrast 
to the previously used techniques, the energy gap is controlled by the film thickness. Our transport 
measurements confirm that the quasiparticle poisoning is suppressed and clear 2e periodicity is 
observed only when the island is made much thinner than the leads. This result is consistent with 
the existing model and provides a simple method to suppress quasiparticle poisoning. 
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Parity effect in a small superconducting island, 
whether it contains an unpaired electron (quasipar- 
ticle, QP) or not, has been extensively studied in 
1990s. 1 - 2 ' 3 ' 4 ' 5 Recently, this topic has attracted research 
interest again 6 ' 7 ' 8,9 because of its importance for the 
quantum bit (qubit) application of the small Josephson 
junction circuits— It is well known that the super- 
conducting charge qubit is least sensitive to the main 
decoherence source, i.e., background charge fluctuation, 
at the charge degeneracy poinb 1 ^ where the electrostatic 
energy of two charge states differing by 2e are equal (e: 
single electron charge). This is because the first deriva- 
tive of the excitation energy of the qubit with respect to 
the reduced gate charge n g — C g V g /2e becomes zero at 
this point n g = 0.5. It should be noted that the curva- 
ture of the energy band, which can be utilized for the 
readoubi 3 ^ or the coupling— of the qubits, is maximum 
also at this point. Thus, it is important to operate the 
qubit at the degeneracy point. As long as the equilib- 
rium states are concerned, operation at the degeneracy 
point can be achieved if E c — Ej/2 is smaller than A, 
where E c , Ej , and A are the single electron charging en- 
ergy, the Josephson energy of a single junction, and the 
superconducting energy gap, respectively (see Fig. |TJi). 
In this case, the ground state is always even-parity (no 
QP) state and its energy is perfectly 2e periodic in n g 
(maximum at half odd integer), which should produce 
2e periodicity of, for example, the gate modulation of 
the supercurrent in the superconducting single electron 
transistor (SET). In reality, however, it is not a trivial 
task at all to observe this 2e periodicity even for the is- 
land having E c —Ej/2 < A probably due to the existence 
of nonequilibrium QPs£ 

Recently, Aumentado et. al. reported 6 " that it is possi- 
ble to reproducibly obtain 2e periodicity of the supercur- 
rent in the superconducting Al SET by making A of the 
island higher than that of the leads (Fig.[IJ)), because it 
works as a barrier, which prevent nonequilibrium QPs in 
the leads from entering the island. They realized such a 
spatial profile of A by introducing a small amount of oxy- 
gen during the evaporation of the island. Similar effect 
has been reported for a Cooper-pair box^ where the re- 
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FIG. 1: (a) Schematic diagram of the energy bands of the 
charge qubit as a function of n g . Dashed lines represent the 
energy bands for the odd parity state, (b) Spatial supercon- 
ducting gap profile, which is favorable for the QP-poisoning- 
free SET device, (c) Scanning electron micrograph of the 
SET sample, whose island and leads are 10 and 50 nm thick, 
respectively. The gate electrode is not seen in this picture. 



quired spatial profile is obtained by a weak magnetic field 
applied parallel to the film. Motivated by these works, 
we investigate an alternative way of creating a proper 
A profile. Our approach is based on the fact that A 
of an Al thin film depends also on the film thickness. 16 
This suggests that we can obtain the spatial profile of 
A by choosing a proper thickness for the island and the 
leads. This technique is simpler than controlling oxygen 
pressure or magnetic field. Moreover, improper choice of 
the thickness may lead to the QP-poisoned devices. Our 
results show that the A profile controlled by the film 
thickness can change the degree of QP poisoning. 

The SET samples were fabricated by electron-beam 
lithography and double-angle evaporation of Al using 
a trilayer resist (methylmethacrylate-methacrylic acid 
(MMA-MAA) /Ge/poly-methyl-methacrylate (PMMA), 
200/20/50 nm thick) on an oxidized Si wafer with pre- 
deposited gold pads. The pressure inside the chamber 
before Al evaporation does not exceed 4 x 10~ 8 Torr, 
and increases up to 1 x 10~ 7 Torr during the evapora- 
tion. No oxygen doping during the Al evaporation was 
performed in this study. The surface of the first Al layer 
was oxidized by 70 mTorr oxygen for 0.7 to 5 min. Fig- 
ure^ shows a scanning electron micrograph of one of the 
SET devices with the thickness of 10 nm and 50 nm for 
the island and the leads, respectively. The bright area is 
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TABLE I: Sample parameters of large Ej/E c ratio samples. 
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FIG. 2: (a) Intensity plot of I s d as functions of V g and V a d- 
The device has 40 nm thick island and 50 nm thick leads, (b) 
Observed superconducting energy gap of Al as a function of 
the film thickness. 



the evaporated Al. The source, drain and gate electrodes 
are connected to the gold pads, which are about 20 
away. No QP traps 5 , are used in the present study. 

In the following, we show two kinds of experiments. In 
the first experiment, we investigate the dependence of A 
on the film thickness by measuring the current-voltage 
(/ — V) characteristic of the SET. For this experiment, 
smaller Ej / E c ratio is desirable for better accuracy in A 
(see below). In the second experiment, we investigate the 
gate modulation of the SET supercurrent with various 
thicknesses. Experimentally, it becomes easier to observe 
supercurrent, if the Ej/E c ratio of the SET is larger. 
Therefore we prepared two sets of SET devices: one with 
small Ej/E c ratio (~ 10 -2 ) and the other with large 
Ej/E c ratio (~ 10 _1 ), by changing the junction area 
and the oxidization time. 

In the first experiment, we use small Ej/E c ratio sam- 
ples. They have typically normal source-drain resistance 
R sd = 260 ktt and E c = 340 fxeV. We fabricate SET de- 
vices with the thickness of the island (ti) ranging from 
10 to 50 nm. The thickness of the leads (ti) is always 50 
nm. Using a 3 He refrigerator at temperatures of ~ 300 
mK, we measure the source-drain current (I s d) as a func- 
tion of the gate voltage (V g ) for different values of fixed 
source-drain bias voltage (V s d), as it is exemplified in an 
intensity plot in Fig. [2i for the sample with U = 40 nm. 
One can see a clear Coulomb blockade diamond struc- 
ture: a sharp rise of the QP tunneling current above the 
Coulomb blockade threshold voltage Vth which is modu- 
lated by V g . From the modulation period, we estimate 
the capacitance between the gate electrode and the island 
C g to be 8.1 aF. The arrays of peaks and dips parallel to 
the diamond structure are due to Josephson-quasiparticle 
cyclesiii As shown in the figure, Vth is modulated be- 
tween 2(A(ii) + A(i;)) and 2(A(U) + A(t t )) + 2E C , where 
A(£) is the superconducting gaps of the Al thin film with 
the thickness t. Thus, we can determine the value of 
A(£;)+A(50 nm) from this measurement^ and A(50 
nm) can be determined from the similar measurement 



Sample 


Rsd [kfi] 


E c [fjieV] 


Ej fcieV] 


periodicity 


10/50a 


52 


145 


27 


2e 


10/50b 


53 


141 


27 


2f:< 


10/30 


38 


139 


37 


2e 


40/50 


32 


114 


40 


2e 


50/10 


47 


160 


30 





for the sample with ti of 50 nm, which turned out to 
be 0.2 meV on average. The derived A(t) is plotted in 
Fig. [2b- It is observed that A increases as t is decreased, 
especially below 20 nm. This is qualitatively consistent 
with the previous reports^ The thinnest sample in this 
study gives 20 to 25% higher A as compared to A (50 
nm). This variation of A is as high as the one achieved 
by oxygen doping in Ref. 0, indicating that we can create 
a similar spatial profile of A in the SET by controlling 
the film thickness. In the next experiment, we check the 
effect of the A profile on the periodicity of the SET su- 
percurrent. 
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FIG. 3: I s d ~ Vsd curves at different gate voltages for (a) 
the sample 10/50 and (b) the sample 50/10. Each curve is 
offset by 0.1 nA for clarity, (c) Gate charge dependence of 
the switching current for all the measured samples. Note that 
random offset charges have not been subtracted. 

Based on the result of the first experiment, we fabricate 
SET devices (large Ej/E c ratio) with several combina- 
tions of ti and ti, by which we will denote the samples 
as U/ti hereafter. Sample parameters are summarized in 
Tabled In all samples, the thinner layer was evaporated 
first. The measurements were performed in a dilution re- 
frigerator at a base temperature of 50 mK using CuNi 
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coaxial cables with RC filters as dc lines. Using the 
current-biased configuration, we swept I s d at a rate of 
~0.1 nA/s and measured V s d for different values of fixed 
V g . Here we show the results of such measurements for 
two samples, which have opposite A profiles: the sample 
10/50a (higher A for the island) in Fig. and the sam- 
ple 50/10 (higher A for the leads) in Fig. 12b. Each curve 
in the panels is offset by 0.1 nA for clarity. In Fig. 
current steps are observed near zero V s d, to which we 
attribute the supercurrent. They are modulated by V g 
and the period is twice larger than that of QP tunneling 
current observed at higher V s( i, indicating that the su- 
percurrent is 2e periodic. In Fig. 03, we observe the QP 
tunneling current modulated by V g , but the supercurrent 
is hardly observed. 

From these measurements, we obtained the switching 
current I s and plotted as a function of the reduced gate 
charge n g = C g V g /2e for all the samples in Fig.^. Here, 
I s is normalized by the theoretical maximum Iq/2, where 
Iq = irA e ff/2eRN is the Ambegaokar-Baratoff critical 
current, Rjy is the normal-state resistance of the single 
junction, and A e// = 2A(t i )A(t I )/(A(t j ) + A(tj)) M- 
We assume that both junctions in the SET have same re- 
sistances. We observe clear 2e periodicity in the samples 
10/50a, 10/50b and 10/30. In all these samples, A of 
the island is about 20% higher than that of the leads. In 
the sample 40/50, we still see the 2e periodicity, but the 
contrast of the gate modulation becomes small. In the 
sample 50/10, the amplitude of the supercurrent peaks 
is very much suppressed and shows no clear modulation. 

The fact that we observe clear 2e periodicity only 
in the samples where A of the island is higher than 
that of the leads is consistent with the previous report. 6 
Namely, in those samples, we create a spatial profile of 



A as shown in Fig. QJ). The energy gap difference \SA\ 
(= \A(ti) — A(tj)|) is about 40 /xeV, which works as a bar- 
rier preventing nonequilibrium QPs created in the leads 
from entering the island. As discussed in Ref.0, observa- 
tion of 2e periodicity does not mean complete absence of 
QPs. In sample 10/50a, 6 A + 5E® 1 becomes positive at 
n g > n g0 = 0.32, where 5Ef = AE c n 2 g - 4E C (0.5 - n g ) 2 
is the electrostatic energy difference between no QP and 
one QP states of the islandiSS It means that QP can 
tunnel into the box in the vicinity of the degeneracy 
point (e.g. 0.32 < n g < 0.68). Although no clear sign of 
poisoning is seen, the small magnitude of 2I s /Iq may be 
the indication of this QP poisoning, 2 -! because our cur- 
rent ramping time is much longer than the inverse of the 
QP tunneling rate, which is reported to be ~ 10 ^tsi^ In 
sample 50/10, one QP state is favorable in a wider range, 
0.13 < n g < 0.87, which leads to the poisoning of the su- 
percurrent. Although we expect the gate modulation of 
I s with le periodicity, it is not clearly observed probably 
because the magnitude of I s is so small that it can be 
easily smeared out by the external noise. 

In summary, by measuring the I — V characteristics of 
the SET, we have shown that A of Al depends on the 
film thickness. We have measured the gate modulation 
of the supercurrent of the SET devices, whose spatial 
profiles of A were controlled by the film thickness. Clear 
2e periodicity is observed only when the island is much 
thinner than the leads, demonstrating that the A profile 
controlled by the film thickness can help to suppress QP 
poisoning. 
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Another reason for the observed small 2I s /Iq is insufficient 
filtering for the external noise in the measurement lines. 



